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We analyze the reactions A + A ~^ III and y4 + _B — ^ occurring in a model of turbulent flow in two 
dimensions. We find the reactant concentrations at long times, using a field-theoretic renormalization 
group analysis. We find a variety of interesting behavior, including, in the presence of potential 
disorder, decay rates faster than that for well-mixed reactions. 
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I. INTRODUCTION 

The behavior of chemical reactions in turbulent flow 
determines how certain types of chemical reactors func- 
tion, how combustion occurs in engines, how smog is 
produced in the atmosphere 0|, and how certain types 
of planktonic predators feed in the ocean |]. Reac- 
tive turbulent flow is usually analyzed by continuum re- 
action/transport equations for the reactants [|l],D. We 
now know, however, that the upper critical dimension in 
which such mean- field equations fail is two . We pro- 
vide here the first renormalization group, field-theoretic 
treatment of reactive turbulent flow in two dimensions. 

We consider a low concentration of reactants immersed 
in an isotropic, turbulent fluid flow. The interplay of 
reaction, turbulent mixing, and trapping by disorder will 
be shown to lead to novel kinetics at long times. So as to 
access the most interesting regime, we will consider a two- 
dimensional system. At low reactant concentrations, the 
dynamics of the fluid will not be affected by the kinetics 
of the reaction. As such, the effect of the turbulence is 
simply to advect and to mix the reactants. 

Our intention is not to derive a theory of turbulence 
but rather to derive a theory of bimolecular kinetics in 
the presence of isotropic turbulence and potential dis- 
order. We, therefore, employ the same statistical the- 
ory of turbulence conventionally used to study turbulent 
transport of passive scalars |^,^. For a review of this 
approach, see |]lO|. We assume, in particular, that the 
turbulent fluid which advects the reactants can be mod- 
eled as a quenched, random, Gaussian velocity field with 
the correct statistics. While the fluid velocity stream- 
lines produced by this conventional approach do not sat- 
isfy the Navier-Stokes equations, the correct transport 
properties of the reactants are captured. The correct 
Kolmogorov energy cascade and Richardson separation 
laws, for example, are produced. One can imagine using 
a more detailed model of the fiuid mechanics. Avellaneda 
and Majda have, for example, used statistical fiow fields 
that depend on both space and time to model turbulent 
transport [pT|-p^. On an even more detailed level, one 
could use statistical flow fields that satisfy the Navier- 
Stokes equations. Renormalization group theories for 
flow fields of this type have been derived by Forster, Nel- 



son, and Stephen and later by Yakhot, Orszag, and 
coworkers [p"5|-pjt. We settle for the simplest descrip- 
tion of the fiuid mechanics that captures the essence of 
turbulent transport. Our results should not be sensitive 
to this assumption, since the physical processes that ap- 
pear to dominate the long-time kinetics depend only on 
the overall transport properties of the fiuid. In support 
of the generality of our results, we note that Deering and 
West have used mean field theory to analyze the reaction 
A -|- i? — > in a time-dependent, but spatially uncorre- 
lated, model of isotropic turbulence. Where mean field 
theory is expected to work (where the renormalized reac- 
tion rate does not enter the prediction for the long-time 
reactant concentration) their results agree with ours. 

In this article, we use a field-theoretic approach to an- 
alyze reaction kinetics in a model of turbulent fluid flow. 
A definition of the fluid flow and the reaction kinetics is 
given in Section II. The reactive turbulence problem is 
mapped to a field theory that is convenient for analysis 
in Section III. The long-time behavior of the A + A ^ ^ 
reaction is derived by a renormalization group analysis 
in Section IV. The long-time behavior of the A + B ^ 
reaction is derived in a similar fashion in Section V. We 
conclude with a comparison to experimental results in 
Section VI. 

II. DEFINITION OF REACTIVE TURBULENT 
FLOW 

In our formulation, an isolated reactant undergoes bi- 
ased Brownian motion in the fluid streamlines, reacting 
at a given rate with other nearby reactants. In the ab- 
sence of reaction, the motion of the particles can be de- 
scribed by a Langevin equation: 

fl-X 

^ = /3Z?F(x,;)+ Vit), (1) 

where the inverse temperature is given by /3 = l/(fcBT), 
and D is the diffusivity. Here the position of particle i, 
Xi, undergoes advection due to forces from the fluid flow 
and diffusion due to forces from the random, thermal 
motion of the fluid. The random, thermal noise has a 
correlation determined by the diffusion coefficient: 
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(2) 



We choose the forces coming from the fluid stream hnes 
so that they mimic tm'bulence. Defining 



(i^^(x)^^,(x')) = G^.(x-x') 



(3) 



we choose 



G^i.(k) = X00(k) {Si_^;,k^ - kf^k^) + x«n(k)fcpfc^ , (4) 

where the Fourier transform of the correlation function is 
G^y(k) = J d''xG^;^(x) exp(tk • x) in d dimensions. Here 
X00(k) is the correlation function of the stream function 
that gives rise to the turbulent fluid flow, and x„tj(k) is 
the correlation function of a quenched, random potential, 
m(x), which we have included for generality. To mimic 
turbulence, we choose 



X««(k) = 



a 



k^+y 

1 

k^+v 



(5) 



Isotropic turbulence is modeled by y = 8/3 and 7 = 0. 

Of course, we are interested in reactive, turbulent flow. 
We consider two different reactions: 



and 



A + A 



A + B 



(6) 



(7) 



Here Aq is the conventional reaction rate. We place our 
reactants on a square lattice, of lattice spacing h. This 
lattice spacing implies a cutoff in Fourier space of A = 
27r//i. Reaction occurs between two particles, at rate 
Xo/h'^, only when they are on the same lattice site. The 
diffusion and advection occurs on this same lattice. 



III. FIELD-THEORETIC REPRESENTATION 

The quantity of interest is the long-time concentra- 
tion of the reactants. The presence of the quenched fluid 
stream lines and the quenched, external potential makes 
direct analysis of the dynamics rather difficult. Pertur- 
bation theory fails due to singularities in the forces at 
small k. We, therefore, map the above description onto a 
field theory and analyze the field theory using renormal- 
ization group theory. We assume that the concentration 
of reactants is initially Poisson, with average density no. 

A field theory is derived by identifying a master equa- 
tion, writing the master equation in terms of creation 
and annihilation operators, and using the coherent state 
representation PJT^]. The random potential is incorpo- 
rated with the replica trick Q, using N replicas of the 
original problem. For reaction (^), the concentration of 



A at time t, averaged over the initial conditions, ca(x, t), 
is given by 



ca(x, i) = lim (a(x,t)) , 



(8) 



where the average is taken with respect to exp(— S*^^), 
Saa = J '^'^^J ^ dtcLai^^t) [dt - DV"^ + Sit)] a„(x,t) 

r-tf 



atidt2 



dt 



2a„(x,t)a2(x,i) 
- no y d'^xaa (x, t) 



'kik2k3k4 

x(27r)'^(5(ki +k2 +k3 +k4) 

xaai(ki, ii)dai(k2, ii)aQ2(k3, i2)aa2 (^4,^2) 
x[ki • (ki -hk2)k3 ■ (ks -hk4)x™(|ki -hk2|) 

-fki X k2 kg X k4X00(|ki -f k2|)] . (9) 



Summation is implied over replica indices. The notation 
stands for / d''k/(27r)''. The upper time limit in the 
action is arbitrary as long as > t. We do not dwell 
on the construction of this field theory. It differs from 
that for reaction in a random potential field only by the 
inclusion of the random stream line terms 

For distinct reactants, reaction (^, a field theory can 
also be derived. The relevant action has the form 

Sab = / d'^x /" ^ dtd„(x,t) [dt ~ DV^ + S{t)] aa(x,t) 



d''x / dt6„(x,t) [at - ZJV^ + (5(i)] 6a(x,t) 



aa(x,t)aa(x, (x,t) 



-f Ao J d'^x j dt 
+bai's.,t)aa (x,t)6Q(x, i) 

+aa (x, t)aa (X, t)ba (x, t)ba (x, t) 

-no J d'^x [da(x,0) -)-6a(x, 0)] 



'kik2k3k4 

x(27r)''5(ki-|-k2 + k3+k4) 
X aQj(ki,ii)dQj(k2,ii) - foa^ (^i, ti)&a2 (^2, ^i) 

003(^3,^2)003(^4,^2) - &a4(k3,t2) 604(^4,^2) 

x[ki • (ki +k2)k3 • (k3 +k4)x™(|ki +k2|) 
-fki X k2 k3 X k4X00(|ki + k2|)] . (10) 

This action also differs from that for reaction in a random 
potential field only by the inclusion of the random stream 
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line terms 0. The concentrations averaged over initial 
conditions are given by 



Cylfx, t) = lim (a(x, t)) 
cb(x, t) = lim (fo(x, <)} 



(11) 



where the average on the right hand side is taken with re- 
spect to exp(— iSab)- So as to reach the most interesting 
scaling limit, we have taken the initial average densities 
to be the same, ca(x, 0) ~ cb(x, 0) = uq. For simplic- 
ity, we have also assumed equal diffusivities of the two 
species. Da = Db = D. Note that the A and B parti- 
cles experience an identical force due to the fluid stream 
lines but an opposite force due to the quenched, external 
potential. 



IV. THE A + A- 



REACTION 



Let us first consider how turbulence will affect the 
A ~\- A ^ % reaction in the absence of potential disor- 
der. Without turbulent mixing, this reaction is diffusion- 
limited. The concentration decays at long times as 
CA{t) ~ \n{t/to)/{^T:Dt), with « h"^ / D |]. This decay 
is slower than the Cyi(i) ^ \/{k*t) that would be pre- 
dicted by simple mean-field kinetics for a well-mixed re- 
action with effective reaction rate k* . Turbulence mixes 
the reactants, tending to eliminate the transport limita- 
tion on the reaction rate. As we will see, turbulence will 
cause the reactant concentration to follow the mean-field 
result, with k* < Aq. 

We analyze the field theory (^) via renormalization 
group theory. The flow equations in two dimensions, to 
one loop order, are 



din tiQ 



= 2 



dl 

din A A 



dl 
din g 



dl 



A-kD 

y-^g, 



(12) 



where the dimensionless coupling constant is given by 
g = aP'^A~y /(An). The dynamical exponent is given by 



z^2-g 



(13) 



We see that the flow equations lead to a non-zero fixed 
point for the coupling g* = y/2. This fixed point is 
probably exact, if we assume that the reaction does not 
affect the transport properties P,p^. 

We determine the long-time decay from the flow equa- 
tions via matching to short-time perturbation theory 
1^. The flow equations are integrated to a time such 
that t{r) = texp[— z{l)dl] = to. At short times, 
we find the mean-square displacement of an unreactive 



particle from (r^(t(r),r)) = ADt{l*) and the concen- 
tration of reactants from CA{t{r*),l*) = l/[l/no(r) -I- 
X{l*)t{l*)]. The long-time asymptotic values are given 



by scaling: (r^(t)) 



(r2(i(r),r)) and CA{t) 



e ^' CA{t{r),r). This procedure gives 



CA{t) 



f 1 



\2-KDy XqJ t 



(7 = 0) 



(14) 



where we recognize the mean field result with effective 
reaction rate 1/fc* = l/{2TrDy) + l/Xo. The mean square 
displacement is given hy{r^{t)) - ADt{t/to)y/^^^y\ 
which is an exact law |10(] in the absence of reaction. 
As expected, turbulence leads to the well-mixed, mean- 
field result for the concentration decay, with an effective 
reaction rate fc* < Aq. We have determined the effective 
reaction rate as an expansion in the parameter y, which 
measures the degree of mixing of the fluid. For simplicity, 
here, and below, we have used used only the fixed point 
value of g when integrating the flow equation for X(l). 

What will happen in the presence of a quenched, ran- 
dom potential? We have seen that fluid streamlines in- 
crease the mixing of the reactants. We have previously 
shown that a quenched, random potential in the absence 
of turbulence leads to a slowing down of the reaction 
P . This occurs because the reaction becomes diffusion- 
limited at long times, and the random potential leads to 
sub-diffusion. We might, therefore, expect that the reac- 
tant concentration for 7 7^ will be lower than that for 
the case of 7 = for arbitrary a. In fact, a subtle trap- 
ping effect due to the potential leads to increased decay 
rates, above that for 7 = 0, for some intermediate values 
of 7. 

The flow equations that result in the presence of the 
random potential are 



din no 
dl 
dlnA 

dl 
d In g^ 

dl 
din ffo 

dl 



X 



AttD 

y - 
= y- 



357 ~ do- 



(15) 



where the dimensionless coupling constants are given by 
g^ = jP^A-y/iAir) and g^ = aP^A-y /(An). The dynam- 
ical exponent is given by 



z ^ 2 + g-^ - ga 



(16) 



These flow equations lead to fixed points for the couplings 
g* = y/2 and g* = (j/ci)y/2. The flow diagram for the 
couplings is shown in Figure |^. Including higher order 
diagrams in the flow equation leads to a bending of the 
fixed line, more strongly for larger 7/(7, with the coupling 
fiows remaining linear [n9| . 
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FIG. 1. The flow diagram for the dimensionless couphngs 
in the A + A —> II) reaction. The flxed hne is shown in bold 
for y — 1. The fixed point reaction rate is finite on the solid 
part of the fixed line and vanishes on the dashed part of the 
flxed line. 



The matching to determine the asymptotic concentra- 
tions leads to two regimes. For weak potential disorder 
with 87 < a, there is no finite fixed point value for X(l), 
and we have 



CAit) 



1 



A7:D{g* ~ 3g; 

^ \ -29;/(2+ff;-g:) 



(37 < <j) 



(17) 



For strong potential disorder with 87 > cr, there is a finite 
fixed point value A* = 47rl?(3(j'* — g%), and we have 



(aX-al)/ {"^+91-91) 



(87 > <y) ■ (18) 



The maximum rate of decay occurs for 87 = rr, in which 
case we have 



CA{t) 



8Tr{l~y/6)Dt 



t-y/id-y)^ (87 = 0-). (19) 



In all cases, the mea-n square displacement is given by 
(r2(i)) ~ ADt{t/to)^3:-g;)/{2-g:+g;)^ 

Examining the asymptotic decay laws ([l7|)-(|l9|), we see 
that a small amount of potential disorder added to the 
turbulent fluid mixing leads to an increased rate of reac- 
tion. As the potential disorder is increased, eventually 
the rate of reaction decreases. The exponent of the con- 
centration decay is shown in Figure ^j. This result is 
rigorously valid for small y and finite values of 7 and 
a. For finite values of y, we expect qualitatively similar 
behavior. 

How can potential disorder, which tends to slow down 
the diffusing reactants, lead to an increased rate of reac- 
tion? The potential disorder creates regions of low en- 
ergy, which tend to attract reactants. The local density of 
reactants in these deep energy wells is significantly higher 




FIG. 2. The decay exponent for the A + A ^ ([} reaction: 
CA{t) ~ (const)f~". The figure is shown for y = 8/3, which 
mimics isotropic turbulence. 



than the average density. The reaction rate in these re- 
gions, therefore, is higher than would be predicted by 
mean-field theory based upon the average density. The 
turbulent mixing flows continuously replenish the reac- 
tants in these wells as the reaction occurs. In this way, a 
reaction rate significantly higher than that for a perfectly 
well-mixed system arises. 



V. THE A + B 



REACTION 



We now turn to the A + B ^ (d reaction. In this case, 
the A and B reactants are attracted to different regions 
of space by the external potential, and so there is no 
mechanism for super-fast reaction. 

The flow equations for this case are 



dlnriQ 

dl 
d\n\ 

dl 
din 

dl 
dlnga 

dl 



A 



47rL> 

y - 

= 2/ - 2g<T 



(37 + 9,y) 



The dynamical exponent is given by 



■9i 



(20) 



(21) 



We first consider the case of no external potential. If, 
in addition, there is no turbulent flow, the A and B reac- 
tants segregate into distinct regions in space. This seg- 
regation leads to a severely diffusion-limited reaction at 
long times. The concentration decays as CA{t) = csi^) ~ 
[no/(87r2Dt)]i/2 Q. Allowing for turbulent mixing, we 
expect the reaction to become more well-mixed, with less 
segregation and faster reaction. In fact, there will be a 
transition to a region that is reaction limited for strong 
enough mixing flows. 



4 



FIG. 3. The decay exponent for the A + B ^ reaction: 
CA(t) ~ (const)i~°'. The figure is shown for 7 = 0. 

We perform the matching to determine the asymp- 
totic decay. In the transport-Hmited regime, we use 
CAit{l*),l*) = [no(/*)/(87r2L)<(r))]i/2. the reaction- 
hmited regime, we use CA{t{l*),l*) = l/[l/no(/*) + 
X{l*)t{l*)l as before. 

For weak fluid mixing, the reaction will be in the 
transport-limited regime, whereas for strong fluid mix- 
ing, the reaction will be in the reaction-limited regime. 
Specifically, for weak mixing we have 



CA 



1/2 



-a/(8-2t/) 



(2/ < 2,7 = 0) . 

(22) 



For strong mixing, we have 
/ 1 



CA{t) 



1 

Ao 



— -, (y> 2,7 = 0) 



(23) 



\2Tr Dy 

The exponent of this decay is shown in Figure |^. Note 
that for isotropic turbulence (jj = 8/3), the reaction is 
always in the reaction- limited regime in two dimensions. 

Potential disorder will slow down the reaction, both be- 
cause reactants are attracted to different regions of space 
and because the transport of reactants to each other is 
slowed. Our one-loop flow equations predict that in the 
transport-limited regime 



CA{t) 



no 



1/2 



(9:-9;)/(4-23;+2g-) 



t 

{g* + 5; < 1) 



(24) 



where, as before, we have 5* = y/2 and g* = (pija^yj^. 
In the reaction-limited regime, we have 



CA(t) 



1 



1 



_47ri?(.g* +.g;) Ao_ 

^ \ 2g;/(2-g;+ff;) 

The exponent of this decay is shown in Figure ^. The 



(l<5:+5:<2) . (25) 



y/CT 

FIG. 4. The decay exponent for the ^4 -|- _B — > reac- 
tion: CAii) ~ (const)t~". The figure is shown for y = 1.2. 
The reaction is transport-limited on the solid curve and reac- 
tion-limited on the dashed curve. The curve is strictly valid 
only for small 7/cr. 



exponent is valid for arbitrary y and small 7/cr. The 
effective reaction rate in Eq. (p5|), however, may contain 
corrections higher order in y. 



VI. CONCLUSION 

Experiments to test our predictions for isotropic tur- 
bulence, Eqs. ( |l^ ) and (|2^), would be relatively simple 
to perform. The behavior of the prefactor would be the 
quantity of interest. Reaction conditions of the type that 
we consider could be realized in reactions between ionic 
species confined to two-dimensional fluid films that are 
surrounded by spatially-addressable electrodes or media 
with ionic disorder that is not equilibrated. The elec- 
trodes or disordered media are necessary to generate a 
potential with the required, singular correlation function. 
The required isotropic turbulence can be generated in the 
standard fashion. Fluid flows less strong than isotropic 
turbulence (y < 8/3) could be observed in regions of de- 
veloping turbulence. 

A recent experiment by Paireau and Tabeling has seen 
an enhancement of the effective reaction rate between 
ions in a chaotically-mixed, two-dimensional, fluid with 
attractors [|2^. In this experiment only the prefactor to 
the reactivity was enhanced. The decay exponent re- 
mained at unity because the disorder was technically ir- 
relevant. We are unaware of experiments, to date, that 
can test our predictions for technically relevant disorder, 
Eqs. dl^-dll) and dD-dl). 
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